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ABSTRACT 

We use spectra of more than 56,000 M dwarfs from the Sloan Digital Sky Survey (SDSS) to create 
a high-latitude extinction map of the local Galaxy. Our technique compares spectra from the stars in 
the SDSS Data Release 7 M dwarf sample in low-extinction lines of sight, as determined by Schlegel, 
Finkbeiner, & Davis, to other SDSS M dwarf spectra in order to derive improved distance estimates 
and accurate line-of-sight extinctions. Unlike most previous studies, which have used a two-color 
method to determine extinction, we fit extinction curves to fluxes across the spectral range from 5700 
to 9200 A for every star in our sample. Our result is an Ay niap that extends from a few tens of pc 
to approximately 2 kpc away from the sun. We also use a similar technique to create a map of Ry 
values within approximately 1 kpc of the sun, and find they are consistent with the widely accepted 
diffuse interstellar medium value of 3.1. Using our extinction data, we derive a dust scale height for 
the local Galaxy of 119 ± 15 parsecs and find evidence for a local dust cavity. 

Subject headings: ISM: dust, extinction — Galaxy: local interstellar matter — stars: late-type — 
stars: low-mass 



1. INTRODUCTION 

The distribution of Galactic dust is highly vari- 
able throughout the Galaxy, and despite many recent 
advances, its precise three-dimensional distribution is 
largely unknown. Because dust both attenuates and red- 
dens sources, the inability to properly correct for it places 
serious restrictions on the precision of many observa- 
tional data products. Dust is also important for star 
formation, since it is thought to be a catalyst for molec- 
ular gas formation and is often found in the cold, dense 
regions where new stars are born. 

Curren tly the most widely used large- s cale d ust study 
is that of iSchlegeL F inkbeiner. fc DavisI ()1998L hereafter 
SFD), which created a two-dimensional map with total 
line-of-sight Milky Way dust co lumn densities calc ulated 
from 100 and 240 /im emission. lArce fc GoodmanI ([l999 ) 
determined that the SFD dust maps overestimate red- 
dening by a factor of 1.3-1.5 in regions of smooth extinc- 
tion with Ay > 0.5 and may underestimate red dening; in 
regions of steep extinction gradients, and Peek fc G ravesI 
f2010) used the reddening of standard galaxies to find 
occasional deviations from SFD of up to 50%. In spite of 
these occasional inaccuracies, the SFD dust map gener- 
ally allows observations outside the Galaxy to be cor- 
rected for the appropriate amount of dust. However, 
sources inside the Galaxy lie in front of only a portion 
of the Milky Way's total dust and thus are likely atten- 
uated and reddened by only a fraction of the total dust 
column. This fraction varies based upon the 3D distri- 
bution of stars and dust in the Galaxy. 

Dust column density can be independently inferred 
from the total line of sight extinction to a given star. 
The extinction affecting stellar spectra can be measured 
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empirically by comparing intrinsic stellar colors to ob- 
served colors, using the formula: 



Av = Rv[E{B-V)], 



(1) 



where Ry is the ratio of total extinction, Ay^ to redden- 
ing, E{B — F), and is an indicator of the dust grain size 
distribution and composition. Ry has an average value 
of 3.1 for the diffuse interstellar medium, and ranges be- 
tween about 2.1 and 5.5 for most sight lines in the Galaxy 
(jCardelH. Clavton. fc M athis 1989, hereafter CCM89). 
Highe r values of Ry are typical of dense regions (e.g. 
Roman-Zuhiga et al.|[2Q07i ). 

The relationship between total extinction and redden- 
ing can be described by an extinction law, usually a 
polynomial fit between the measured color excesses (e.g. 
E{B — V)) of different bands normalized to Ay. The 
most commonly used extinction law is that of CCM89, 
which fits a law of the form A{X)/Ay = a{x)^h{x) / Ry to 
extinction measurements in Johnson photometric bands, 
where a{x) and h{x) are polynomial functions of x, which 
is in inverse wavelength units (x = 1/A fim~^). In 
CCM89 and other similar extinction laws (e.g. Rieke 
& Lebofsky 1985; O'Donneh 1994; Fitzpatrick 1999), 
the shapes of extinction curves for optical and longer 
wavelengths depend only on the amount of dust and 
Ry^ but their i?y-dependence varies considerably with 
wavelength. Infrared curves depend only marginally on 
Ry while shorter- wavelength optical curves have a strong 
i?y-dependence. 

Most of the extinction laws discussed above determine 
extinction by assuming that the reddening over each 
Johnson photometric band corresponds to reddening at 
the effective wavelength of the filter. However, because 
reddening increases toward shorter wavelengths over a fil- 
ter's width, the effective reddening wavelength is slightly 
shorter than the filter's effective wavelength. By us- 
ing modified effective reddening wavelength s as the basis 
for a spline fit to extinction values, iFitzpatrickl () 19991 ) 
presents a reddening law that is more accurate when 
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used with narrow-band photometry, other photometric 
systems, and spectroscopic data. 

In order to overcome the shortcomings of a two- 
dimensional extinction map (e.g. Burstein & Heiles 1978, 
SFD), several attempts have been made to use color 
excesses and extinction laws to m odel Galactic e xtinc - 
tion and dust in three dimensions. iHakkila et al.l ()1997i) 
completed one of the most comprehensive all-sky empir- 
ical extinction models to date by combining; the extinc- 
tion surveys of Fitzserald' fl968'),'Necke]jLKlare' ('1980'), 
[Berdnikov & Pavlovskaya (1991), and Arenou et al] 
([l992f ). Using Equation 1 (choosing an average Ry of 
^3.1) in conjunction with intrinsic and obs erved absolute 
B and V magnitudes. IHakkila et all ()1997l ) computed the 
line of sight extinctions for around 78,000 stars. 

These early models were important for providing an 
initial estimate of 3D extinction, but suffered from lim- 
ited resolution and large measurement er ror due to un- 
known Rv values. The combined model of lHakkila et al.l 
(|T997[ ) was effective at identifying absorbing clouds 
within 5 kpc, but was only intended to identify large- 
scale extinction near the Galactic plane. The stud- 
ies that comprise the comb ined model w er e fai rly lim- 
ited, and even though the I Arenou et al.l ()l992f ) study 
had the greatest number of stars, its accuracy was lim- 
ited to about 1 kpc wi th a resolution o f sever al de- 
grees^ jFitzgerald" (T968^, iNeckel fc Klard ([19801 ). and 
iBerdniko v & Pavlovskaya (j 199 If ) also had limited reso- 
lution due to the rare hot (mostly B-type) stars that 
comprised most of their samples. There were also signif- 
icant differences in data quality between the studies in 
the combined model, leading to poorly defined errors in 
many areas of the map. 

The assumption that Ry equals 3.1 introduced further 
error in Hakkila et al. (1997) due to the short opti- 
cal wavelengths used to compute extinction. According 
to SFD, extreme Ry values can cause deviations from 
the Rv = 3.1 extinction law of up to 30% in the V 
band and 55% in the B band, so E{B — V) could be 
subject to considerable error. Although Ry appears to 
be approxima t ely 3.1 for most sight li nes in the G alaxy, 
ISchlaflv et all ^B) and Schlafl v Finkbeineii (|2010 ) 
found a broad distribution of Galactic Ry values. 

In recent y ears there have been considerable improve- 
ments to the IHakkila et al.l (jT997) map at low Galactic 
latitudes, espec ially looking towards the Galactic center. 
[Marshall et al.l (p306) created an extinction map using 
data from the Two Micron Ah Sky Survey (2MASS) for 
b < 10° and I < 100°, an area that is often ignored 
in extinction maps. The iMar shall et al.l ()2006f ) model 
of the inner Galaxy contains high-resolution extinction 
values for a portion of the sky close to the Galactic 
plane. They compared J — Ks colors to predicted in- 
trinsic colors from , the Besancon Model of the Galaxy 
(jRobin et al.l l2003l ) with K and M giants to determine 
the extinction in the inner Galaxy. Because infrared ex- 
tinction curves are assumed to be independent of 
this approach is more robust to Ry variati on than those 
using B — V color excesses. Similarly, Sale et al l (|2009) 
created an algorithm capable of determining extinction 
in A to early K-type stars by comparing intrinsic to ob- 
served r—i colors (wavelengths with low i?y-dependence) 
from the INT/WFC Photometric Ha Survey (^Drew et al.l 
[2005,), which they used to map extinction in high reso- 



lution at low Galactic latitudes. Using a similar method 
as Marshah et al. (2006|, also with 2MASS photometry, 
Gontc harovi (j2009i, .20100 used 0-F type stars to map ex- 
tinction within 1600 pc of the sun at a range of Galactic 
latitudes with greater resolution than is afforded by the 
more rare K and M giant stars. Although higher than 
the resolution of iMarshall et al.l (|2006f ) in local regions, 
the resolution of this map is currently limited to Galac- 
tocentric cartesian boxes of 100 pc per side. 

Photometric and spectroscopic data from th e Sloan 
Digit al Sky Survey (York et al. 2000; Abazaj ian et al.l 
I2009L SDSS) has recently been used to model red- 
dening at higher latit u des. ISchlaflv et al] (|2010f ) and 
ISchlaflv fc Finkbeinerl ()2010f ). with photometry and 
spectroscopy, respectively, used the colors of main- 
sequence turn off stars to test reddening from the dust 
maps of SFD and to provide the first large-scale, two- 
dimensional map of Ry in the Galaxy. Because this sur- 
vey needed to use bright stars with relatively invariant 
photometric colors, the stars used to calculate extinction 
were behind the entire dust column except at latitudes 
near the Galactic plane. 

Missing from these studies is a comprehensive, local, 
three-dimensional extinction map of higher Galactic lat- 
itudes. M dwarfs, as the most ubiquitous stellar types in 
the Galaxy, are the most suitable candidates to pursue 
such a study. However, because of their low luminosity, 
M dwarfs have been almost completely ignored by pre- 
vious extinction studies. In addition, fitting extinction 
curves to photometric data for these stellar types (and 
many others) is much less robust than a spectroscopic 
method due to the fewer data points over which fitting 
can occur. However, with the advent of large spectro- 
scopic surveys such as the SDSS, there are now spectro- 
scopic samples th at contain more than 70,000 M dwarfs 
(jWest et al.l 1201 ll ). Using these data, it is now possible 
to fit extinction curves to a large spectral range. We 
now have the ability to pursue a more precise three- 
dimensional extinction map for high Galactic latitudes 
within about 2 kpc of the sun using M dwarfs. 

Using SDSS data to more precisely model extinction 
also gives us the ability to derive the three-dimensional 
distribution of dust in the Galaxy. The main structural 
features of dust in the local Galaxy are the local bubble 
and an exponential disk. The local bubble has been mod- 
eled in detail by measuring the lack of neutral ga s within 
^50 t o 200 pc of the sun. Most recently, Lallemen t et al.l 
(|2003f ) modelled the local bubble in three dimensions by 
using absorption maps of Nal. However, although the 
neutral gas profile of the local bubble is well-constrained, 
the dust column densities are not. A local underdensity 
of reddening has been established (e.g. Knude & H0gi 
1999) and some of its contours have been shown to cor- 
relate with the gas cavi ty using reddening fr om stars with 
Hipparchos parallaxes (jVergelv et al.ll2QlQ[ ). but the pre- 
cise shape and size of the dust cavity is not well known. 

The scale height of dust in the local Galaxy is also 
relatively unconstrained. The theoretical model of 
Dri mmel fc Spergel (2001) finds a value ~185 pc near 
the sun and the HI model (l ong thought to be a tracer of 
dust) of iKalberla fc KerpI ()2009i) finds an approximate 
value of 150 pc at all Galactocentric radii, but these 
studies have a high uncertainty in the local environment. 
Because M dwarfs are the most numerous population of 
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stars and have be en catalogued in ^eat abundance in 
the nearby Galaxy (jBochanski et al.ll2QlQf ) , they are ideal 
candidates to study the dust abundances in greater detail 
within several hundred pc of the sun. 

In this study, we used an SDSS spectroscopic sample 
of more than 56,000 M dwarfs to model Galactic dust 
and extinction at high latitudes in the local Galaxy. In 
section 2 we discuss our spectroscopic sample, and in 
section 3 we discuss creating our extinction, and 
dust models. Our results are presented in section 4 and 
discussed in section 5. 

2. DATA 

SDSS data have proven to be incredibly valuable for 
large-scale Galactic surveys due to their high quality 
photometry over large solid angles. Ivezic et aTl (j2008 f) 
used SDSS spectroscopy of 60,000 F and G-type main- 
sequence stars to model t empe rature and metallicity in 
the Galaxy. Juric et al.l ()2008f ) used photometric data 
from ~48 million stars to create an unprecedented model 
of the stellar number density distribution. SDSS has also 
been an id eal tool for stu dying the properties of low- 
mass stars. iBochanski et al. (2010) used around 15 mil- 
lion low-mass dwarfs to measure the luminosity function 
and mass functions for these stars. It also enabled large- 
scale studies of magn etic and chromospheric activity 
(jWest et al.ll2004ll2008f ). and the use of low-mass stars as 
kinematic tracers of Galactic d ynamics (Bochanski et al. 
I2007ai [20TTI : iFuchs et al.l[200l Pineda et al. in prepara- 
tion). 

Our spectroscopic sa mple was selected from , the SDSS 
Data Release 7 (DR7; 'Abazaiian et all l2009f ) M dwarf 
spectroscopic catalog (West et al. 20lT|). SDSS spec- 
troscopy was carried out by twin fiber-fed spectrographs 
collecting 640 simultaneous observations. Typical ex- 
posure times were '^15-20 minutes, but exposures were 
subsequently co-added for exposure times of ~45 min- 
utes, producing mediu m resolution spectra with R ~ 
2000 (|York et al.ll2Q00f ). SDSS spectroscopic plates each 
contained 16 spectrophotometric standard stars, which 
were selected by color to be F subdwarf stars. The SDSS 
spectroscopic fluxes were calibrated by comparing the 
standard stars to a grid o f theoretical s pectra from Ku- 
rucz model atmospheres (iKuruczlll992 f) and solving for 
a spectrophotometric solution on each plate. The SDSS 
spectrophotometry matches the PSF photometry to 4% 
rms (Adelman-McCarthy et aP l2008l ). For more infor- 
mation on SDSS spectrophotometric calibration, see the 
discussi ons in the SDSS second and sixth data release 
pape rs (jAbazajian et al.ll2004HAdelman-McCarthv et al.l 
[2OQ8D. 

Our initial sample consisted of 70,823 M Dwarfs (MO 
to M9), which were chosen based on typical M dwarf col- 
ors {r - i> 0.53 diud i - z> 0.3; West et al. 2008 ) and 
had spectral types verified by eye. To ensure the high- 
est quality spectroscopic data, stars with signal-to-noise 
ratios near Ha of less than 4 were removed. Stars with 
colors indicating a possible M dwarf- White dwarf binary 
were also removed (Smolc ic et al. 2004), as were stars 
that didn' t pass th e SDSS photometric quality fiags de- 
scribed in IBochanski et al.1 ()2010f ). The resulting sample 
contained 55,719 stars. 

The radial velocities were measured for all stars in 
the SDSS by cross-correlating high signal-to-noise ratio 



templa tes with each M dwarf spectrum (jBochanski et al.l 
l2007bD . Distances were measured for all stars us- 
ing the photon i etric parall ax technique des cribed in 
IBochanski et~ZI (|20Toh . See lWest et J] (|20TTh for more 
details on the sample selection and the value added quan- 
tities measured for the SDSS DR7 M dwarf sample. 

3. ANALYSIS 

In order to measure the alteration of stellar spectra by 
dust, we selected a set of template spectra from regions 
of the sky where the total integrated extinction (as de- 
termined by SFD) was less than 0.03 magnitudes in the r 
band. We then separated the total set of template spec- 
tra by spectral type. For spectral types M7 and earlier, 
we also put spectra into 5 bins based on the strength of 
the Ti05 molecular band head. The TiO bins were cre- 
ated to minimize the effects of metallicity on extinction 
curve fits (changes in metallici ty can reduce the depth 
of th e TiO molecular features; iGizisI IT997I : iLepine et al.l 
[2007h . The bins were computed so that each one had 
roughly the same number of template spectra, allow- 
ing every star to be compared against several templates. 
Due to the weaker Ti05 features of M8 and M9 stars 
(caused by the condensation of TiO i n their cooler at- 
mospheres; lAllard fc Hauschildtlll995l ) and the scarcity 
of low-extinction template stars for these types, we did 
not separate M8 and M9 templates into Ti05 bins. 

We visually inspected all of the template spectra for 
missing data or low signal-to-noise ratios over parts of 
their spectrum, finding that many of the stars in the orig- 
inal template sample were not suitable for use as "ideal" 
M dwarf spectra. However, we managed to identify ~5- 
10 extremely high quality spectra (out of ~60 originally) 
for each Ti05 bin. For the templates that we selected 
by eye, the median extinction given by SFD was 0.028 
magnitudes in the r band with a median distance of 141 
pc. Each of the template spectra were corrected for their 
radial velocity and spline-fitted to a wavelength array 
that was spaced 69.1 km s~^ apart in velocity (SDSS 
wavelength spacing). 

Due to the low signal-to-noise ratios in the blue por- 
tion of many M dwarf spectra (caused by their red colors 
and low luminosities) and the high i?y-dependence of 
shorter- wavelength extinction laws, we limited our ex- 
tinction fits to the region spanning from 5700 to 9200 
A. For our initial model we assumed Ry was equal to 
3.1 everywhere, which is a reasonable approximation for 
most lines of sight (CCM89). Our assumed Ry does 
not greatly increase our uncertainty because extinction 
at longer wavelengths, especially those longer than the 
V band (~5500 A), has a relatively low i?y-dependence. 
We derived a relationship between fiux, distance, and ex- 
tinction at each wavelength by using the fact that the ra- 
tio of fiuxes between the program star and template star 
is a product of both the distance modulus and the extinc- 
tion. Using the i?y=3.1 extinction curve of Fitzpatrick| 
(1 19991), we derived the following equation, which relates 
the ratio of fiuxes between each star and low-extinction 
template star to the distance and total V band extinc- 
tion: 



2.5log 



/(A) 



/t(A) 



= 5log 



-4(A), 



(2) 



where / is the fiux of the program spectrum, ft is the 
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Fig. 1. — Top: An MO spectrum with some reddening. Middle: 
An MO template spectrum with less than 0.03 magnitudes of r band 
reddening. Bottom: 2.5 times the logarithm of the flux ratios of 
these two spectra, which is equal to the extinction added to the dis- 
tance modulus. This ratio has been flt with the Fitzpatrick ( 1999) 
Rv =3.1 extinction law (green line). The reddening is apparent in 
the slope of the extinction law flt, which had an Ay of 0.88 magni- 
tudes. The SDSS flux is in units of 1 x IQ-^'^erg cm-'^ A-^. 

flux of the template spectrum, d is the distance to the 
program star, dt is the distance to the template star, 
and A(X) i s the extinction at wavelength A given by 
iFitzpatrickl ^999). 

For each program M dwarf in the IWest et al.l (j2Qlir ) 
sample, we corrected for radial velocity and interpolated 
the flux to the wavelength array of our template spec- 
tra. Each program star flux array was divided by the 
flux array of every template star with the same spectral 
type and Ti05 strength. This gav e us values for the 
flux ratio at every wavelength. The IFitzpatrickl () 19991 ) 
extinction curves were parameterized by the IDL pro- 
gram fm.unrecO, which we modified slightly to return 
the A{X)/ E{B — V) extinction curve. Dividing this curve 
by Rv = Av/E{B — V) gave us the relation normalized 
to A{X)/Av. Multiplying the extinction curve by Ay 
returned A(A), which can be described by the equation: 



AiX) = {Av/Rv)[A{\)lE{B - V)] 

Fitzpatrick -> (3) 

where Ay was the only unknown. 

As suming that the derived distances (jBochanski et al.l 
I2Q1QI : West et al. 2011) were correct for the case of low- 
extinction templates, we were left with 2 free param- 
eters, the program star distance, (i, and the V band 
extinction. Ay. We then performed a two-parameter 
Levenber g-Marqu art least-squares fit using mpf itfun in 
IDL (Ma rkwardtl 12009). which returned the best fit Ay 
and distance based on the fluxes at every wavelength 
(Figure [T]). Any residuals in the fits are likely due to 
noisy spectra or an imperfect match of program star to 
template star. Taking the median of the fits for each 
template star gave the stellar distance and absolute V 
band extinction for each program star. The dispersion 
in the values measured for each template star was used 
for the uncertainty. The fits were constrained to have 
both distance and Ay > 0.0 and were weighted by the 
noise at each wavelength. We were able to fit Ay to a 
median uncertainty of 0.043 magnitudes and distance to 

^ Part of the IDL Astronomy User's Library at: 
|htt p : / / idlast ro . gsf c . nasa. gov / 1 



a median uncertainty of less than 5%. This adds to the 
uncert ainty in the original distances given bv lWest et al.l 
(|2Q11[ ). which is typically ~20%. With the new distances, 
we derived a set of cylindrical Galactocentric coordinates 
for each star in the sa mple and hav e updated the sample 
distances in the West e t al.l ()2011l ) catalog. To compute 
the error in Ay, we took the dispersion in Ay for each 
program star and divided by the square root of the num- 
ber of template stars against which each program star 
was fit. 

3.1. Fitting Ry 

A subset of our sample was used to examine the dis- 
tribution of Ry. For the program and template stars 
with signal-to-noise ratios near of greater than 4, we 
used the same procedure described above to perform a 
3-parameter fit, this time including Ry in Equation 3 as 
a free parameter. Due to the high signal-to-noise ratio 
at the blue end of the spectrum for these stars, we were 
able to fit over the larger wavelength range of 4300 to 
9200 A for 9102 M dwarfs (Figure [2]). As with our two- 
parameter fit, the residuals are likely due to an imperfect 
match of program spectrum to template spectrum. A 
plot of the Figure [2] extinction curve residuals is shown 
in Figure [3l The residuals shown are generally small, 
having a median value of 0.037, but tend to increase to- 
ward bluer wavelengths. Unlike with the 2-parameter 
fits described above, where Ay was restricted to have a 
minimum value of 0, we found that mpf itfun was better 
able to constrain Ry when we allowed Ay to be slightly 
negative. We therefore allowed Ay to be as low as -0.036 
magnitudes, the systematic error from using a template 
star with 0.03 magnitudes of r band extinction (assuming 
Ry of 3.1), which allowed robust fits of Ry that would 
be less affected by the systematic error in Ay. 

We found that Ry is a sensitive parameter to the un- 
certainty in the spectrum, due to the fact that it de- 
scribes the subtle changes in the curvature of extinction 
laws. Many of the noisier fits returned values that were 
too high or too low to be plausible for Galactic sightlines. 
We expect that these values, rather than being indica- 
tive of unlikely Galactic dust properties, could be due to 
the effect of random noise or possibly molecular features 
which could be slightly different in the program star than 
in the template star. Our Ry map also demonstrates 
that stars at greater distances, which are more likely to 
have lower-quality spectra, show the greatest dispersions 
in Ry, indicating that these values might not be physi- 
cally meaningful. In addition, the full width at half max- 
imum of our distribution is ~2, so these values appear to 
be well outside the normal range. For these reasons, we 
constrained our analysis to Ry values between 2.1 and 
5.5 (CCM89). Values in this range were returned for 6338 
of the M dwarfs (~70%) when compared with at least one 
template star. The remaining 30% of stars reported un- 
likely Ry values when compared with all template stars. 
The values between 2.1 and 5.5 were kept and used for a 
Galactic model of Ry and to glean a rough sense of the 
overall distribution of Ry values. 

As with Ay, we computed the error in Ry values by 
first taking the dispersion of the median for the Ry values 
reported for each program star. We then divided by the 
square root of the number of template stars against which 
each program star was fit. However, we often found that 
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Fig. 2. — Top: a M3 spectrum with some reddening. Middle: a 
M3 template spectrum with less than 0.03 magnitudes of r band 
reddening. Bottom: 2.5 times the logarithm of the flux ratios of 
these two spectra, which is equal to the extinction added t o the 
distance modulus. This ratio has been flt with the Fitz patricld 
(HMD extinction law (green hue; Ay = 0.81; Ry = 3.31). The 
SDSS flux is in units of 1 x IQ-^^erg 5"^ cm~2 
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Fig. 3. — A plot of the residuals for the extinction curve flt in 
Figure O The residuals are generally small, with a median value 
0.037, but tend to increase towards bluer wavelengths. 

only one template star for a given program star returned 
a realistic value of Ry- For these stars, we used the 
formal 1-a errors reported by mpf itfun. We were able 
to obtain a median uncertainty in Ry of 0.42. 

3.2. Fitting to an Exponential Model 

In an effort to discover the large-scale distribution of 
local dust, we subsequently fit our extinction values to an 
integrated exponential dust model scaling in Galactocen- 
tric radius {R) and height (Z) and uniform with respect 
to the Galactocentric cylindrical angle 0. By doing this, 
we derived the approximate scale height and length of 
the dust. Extinction is related to dust column density 
by the formula: 



74y(x) = / pd{s)ds, 

Jo 



(4) 



whe re is the mean opacity of the dust (jPrimmel et al.l 
|2003). The mean opacity is estimated to be 8.55 x 10 
cm^g~^ by the model of IWeingartner fc Draind ()2001l ) 
with slightly reduced abundances to match the canoni- 
cal Ay/Nn relation and optical constraints from iDraind 
(j2003a .b.). The scale height derived from this model, 



however, does not depend on which value of tzy we 
use. We assumed an exponential dust model scaling 
in the Galactocentric R and Z coordinates of the form 
S(R,Z) = Soexp(-i?/i^i?)exp(-Z/i^z), which we in- 
tegrated analytically to find the extinctions predicted by 
our model. In order to derive the scale height of the 
dust, the integrated extinc tion was fit to th e observed 
extinction by mpf it2Dfuii (jMarkwardtl 120091 ) with Hz, 
Hr^ and So as free param eters. The corrected photo- 
metric parallax distances (Bochansk i et al.ll2010f ) used in 
this fit had a typical uncertainty of ~25%. 

3.3. Fitting to the Local Bubble 

Inspection of our Ay model revealed a potential ab- 
sence of dust near the Galactic plane that may be asso- 
ciated with the local bubble. Although the local bubble 
has only been well-constrained by its absence of neutral 
^as, neutral gas tends to be a tracer for dust (jPraind 
l200 3a) and thus we expect the local bubble to contain 
a significant underdensity of dust. We therefore decided 
to add the local bubble to our exponential dust model. 
Although the shape of the local bubble is highly irregu- 
lar, our limited resolution at small distances is only suf- 
ficient to roughly constrain the size of the bubble (due 
to the small number of nearby M dwarfs that don't satu- 
rate the SDSS detectors). Because of this constraint and 
the approximate shape of the three-dimensional model of 
iLallement et al.l ()2003l ). we decided to adopt a spherical 
model for the local bubble. This model roughly matches 
the shape of the local bubble of ionized gas. 

Using the assumption of a spherical local bubble, 
we adopted an iterative approach wherein we used 
mpf it2Dfun (M arkwar dt 2009) to fit to the scale height 
and scale length of an exponential model with a spherical 
cavity around the sun of varying radii. Our free param- 
eters were the scale height, the scale length, the central 
density (Sq), and the fraction of dust in the local bub- 
ble relative to the fraction predicted by our exponential 
disk model. We first computed the reduced chi-square 
for bubble radii with the dust fraction parameter fixed 
at 0.0. We adopted the bubble radius with the lowest 
chi-square value. Next, we included the dust fraction as 
a free parameter in order to probe the density of dust 
within the local bubble. 

4. RESULTS 

Using the methods described above, we measured Ay 
and distance using 55,719 M dwarf spectra. For more 
than 6,000 high signal-to- noise spectra, we also estimated 
Ry. Histograms of our derived values for Ay and Ry are 
shown in Figures H] and [5l respectively. Using these data, 
we created maps of Ay and Ry^ as well as a best fit model 
for the Galactic dust distribution. We also compared our 
values to those of SFD for all lines of sight. The details 
of our results are discussed below. 

4.1. Extinction Map 

Figure [6] shows the Ay map from our analysis, created 
using the median Ay values in 25 pc by 25 pc Galac- 
tocentric bins compressed in the ^-dimension. As these 
bins are in Galactocentric coordinates, their angular size 
varies with distance while their spatial size remains con- 
stant. The extinction map extends to approximately 2 
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Fig. 4. — A histogram of our values for Ay. The values where 
Av = 0.0 comprise ~40% of our sample (~25,000 stars). The 
distribution falls off exponentially toward higher Ay^ which we 
would expect given that the histogram is dependent on the dust 
distribution, which we expect to scale exponentially with distance 
(from our analytical integration of Equation 4), and the sample 
density of M dwarfs as a function of distance. The sampled M 
dwarf density begins to drop off exponentially after ~200 pc, well 
before the median sampled M dwarf distance of ~500 pc. 



400 




Fig. 5. — A histogram of our Ry values. The distribution peaks 
close to the widely accepted diffuse interstellar medium value of 
3.1, although our median Ry value of 3.38 is located slightly to 
the right of the peak. Regions of higher dust column density are 
expected to account for the distribution's elongated tail. 

kpc above and below the Galactic plane. In the northern 
hemisphere, the distribution of extinction across different 
lines of sight appears to be fairly regular, whereas in the 
southern hemisphere it appears patchier due to the rela- 
tively few lines of sight in the southern SDSS survey area 
(which likely does not represent the average line of sight 
extinctions in the southern Galactic hemisphere). Al- 
though extinction generally increases at larger distances, 
there may be a selection effect at distances greater than 
around 1 kpc from the sun, where proportionally more of 
the stars observed are from low-extinction lines of sight. 
This bias might arise because stars in high-extinction 
lines of sight are generally noisier and less likely to re- 
main in our sample after our signal-to- noise cuts. 

The median derived extinction for the entire sample 
is 0.037±0.029 magnitudes. A histogram of our Ay val- 
ues (Figure H]) shows that the number of extinction values 
measured tends to decrease exponentially towards higher 
Ay. This distribution is anticipated due to the fact that 
Ay depends on the dust distribution, which scales expo- 
nentially along a line of sight, and the M dwarf sample 



density, which also scales exponentially in distance (out- 
side of ~200 pc). This histogram also shows that 43% 
of our sample reported an Ay of 0.0. Out of these val- 
ues, 24% reported no dispersion in their values due to 
the vast majority of the template star fits all reporting 
an Ay of 0.0. SFD report a median Ay of 0.094 mag- 
nitudes for these stars, while they report a median Ay 
of 0.097 for the entire M dwarf sample. The stars to 
which we measure an Ay of 0.0 have a median distance 
of 470 pc, well beyond the typical scale scale height of 
the dust disk, and so many of these values could partly 
be caused by random and systematic error or unusual 
lines of sight. Out of the stars with a detecteable Ay^ 
the median dispersion is 0.043 magnitudes. We assumed 
that these dispersions represent the random errors in our 
extinction fitting; every program star was fit against mul- 
tiple templates and the uncertainly can be gleaned from 
the deviation of all of the fits for a single program star 
from the median value. 

In addition to the random error in our sample, there is 
also some systematic error that arises from the fact that 
our templates are not from perfectly extinction-free lines- 
of-sight. Our template stars were chosen from regions of 
total column extinction of up to 0.03 magnitudes in the 
r band. Assuming an Ry of 3.1, this translates to 0.036 
magnitudes in Ay . For the templates that we selected by 
eye however, the median total column extinction as given 
by SFD was 0.032 magnitudes in Ay, but with a median 
distance of only 141 pc. Our results indicate that this 
median distance is only slightly above the average dust 
scale height of the Galaxy (see Section 4.3), and thus we 
expect that the actual Ay to our templates is less than 
the maximum value. The small median distances are 
due to the template spectra being selected to have high 
signal-to-noise ratios. This process preferentially selects 
brighter stars, which are more likely to be at smaller dis- 
tances and have Ay values lower than the total line of 
sight values. In addition, the difference in intrinsic tem- 
plate spectra extinction could contribute to the reported 
error in our fitting procedure. Because template spectra 
have varying amounts of intrinsic Ay^ the deviation in 
our template star fits will be in part due to the system- 
atic error. Because of this, our systematic error is likely 
less than the maximum value of 0.036 magnitudes. 

Figure [71 shows a map of the error in our Ay map 
(Figure [6]). For the 25 pc by 25 pc Galactocentric bins 
in this map, we computed errors by taking the deviation 
from the median value of Ay for all the stars in each bin 
and dividing by the square root of the number of stars. 
For bins with only one star, we used the error measured 
by our fitting procedure. As expected, the error in our 
values tends to increase at greater distances as spectra 
become noisier and less reliable. For the 25 pc by 25 pc 
bins, the median error is 0.046. The small increase in the 
median error when using these bins likely demonstrates 
the patchiness of the Galactic dust distribution (we have 
compressed all lines of sight in the ^-dimension). 

In addition to the map of Ay in Galactocentric cylin- 
drical coordinates, we also present a map of our sample 
in Galactic latitude and longitude. Figure [8] shows a map 
of Ay that includes all the stars in our sample averaged 
over bins of 4 square degrees. This map shows that the 
SDSS data in the southern hemisphere come from only 
a few lines of sight, which explains why the extinction 
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Fig. 6. — Median Ay (color coding) in Galactocentric cylindrical coordinate (R, Z) bins of 25 pc by 25 pc. The map is flattened in the 
'-dimension. 
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Fig. 7. — Median error (color coding) of our Ay map in Galac- 
tocentric cylindrical coordinate (R, Z) bins of 25 pc by 25 pc and 
flattened in the ^-dimension. Generally, the error in our values 
increases at greater distances. 

in the southern hemisphere of our Galactocentric map 
(Figure [6j) is much more patchy than the extinction in 
the northern hemisphere. 

4.2. Rv Map 

Using the 6,338 stars for which we measured reahstic 
Rv values, we created a map of our median Ry for each 
50 pc by 50 pc Galactocentric cylindrical coordinate bin 
(Figure [9]). As with our Ay map, the bins have an angu- 
lar size which depends upon the distance at which they 
are located. Accurate Ry values proved challenging to 
determine at large distances. Therefore our map of Ry 
only extends to 1 kpc from the sun. When possible. 




Fig. 8. — A Galactic latitude and longitude map of our Ay values 
(color coding) for all the stars in our sample averaged over bins of 
4 square degrees. 

we measured the errors in Ry by the same procedure 
described above for extinction values. When only one 
template star returned a plausible value for Ry, we used 
the fo rmal 1-a errors returned by mpf itf un (jMarkwardtl 
l2009l ). Figure [5] shows the distribution of Ry^ which 
peaks at a value of 3.1, in agreement with the gener- 
ally accepted diffuse interstellar medium value (CCM89). 
Our median Ry is 3.38 with a median uncertainty of 0.42. 
Figure [To] shows Ry plotted as a function of Galactic lat- 
itude. There appears to be an increase in Ry values at 
higher Galactic latitudes, which is contrary to the ex- 
pectation that Ry is higher in the dense r dust of the 
Galactic plane (jRoman-Zuhiga et al.ll2QQ"7l ). 
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Fig. 9. — Map of Ry values (color coding) in Galactocentric 
cylindrical coordinate (R, Z) bins of 50 pc by 50 pc and flattened 
in the ^-dimension. The data give realistic values for Ry to a 
distance of about 1 kpc. 

We also examined the systematic error in our sam- 
ple due to the fact that our template spectra have an 
unknown amount of reddening (up to 0.036 magnitudes 
in Ay)' Because of this, their intrinsic reddening con- 
tributes to the derived value of Ry- We examined this 
error by using the IDL program fm.unred to artificially 
redden our template spectra by an additional 0.036 mag- 
nitudes in Ay using the .Fitzpatrick (199 9) parameteri- 
zation. Our distribution, including the peak at an Ry of 
3.1 was mostly unchanged, but the median value shifted 
from 3.38 to 3.35 when the spectra were reddened. Al- 
though we expect that our template spectra have less 
reddening than the maximum for the reasons discussed 
in section 4.1, there is still some systematic error. There- 
fore, due to the unknown amount of intrinsic reddening 
of our template spectra, the true median Ry could be 
slightly higher than the value we measured. 

We next examined whether the observed distribution 
of Ry was consistent with a Gaussian with the width 
equal to our uncertainty in Ry. To do this, we cre- 
ated a random Gaussian distribution with 9,102 values 
of Ry (our initial sample size) centered at 3.1 and with 
a = 0.42 (our median error). We removed all values that 
were less than 2.1 or greater than 5.5. We performed 
a Kolmogorov- Smirnoff test on the observed and simu- 
lated Ry distributions, and found a 0% probability that 
the two datasets were drawn from the same parent sam- 
ple. The reported distribution of Ry is therefore broader 
than can be explained by a single Gaussian distribution 
with (J = 0.42, and thus is indicative of a real spread in 
Galactic Ry. However, some of the spread is likely due 
to our measurement uncertainties. 

4.3. Fitting to a Simple Dust Model 

To fit an exponential dust model with a local bubble to 
our extinction results, we first reduced the measured un- 
certainties by putting the stars into Galactocentric cylin- 
drical (R, Z, Rsin^) three-dimensional boxes of 25 pc on 
a side at all distances. We weighted each bin by the num- 
ber of stars it contained and only used the stars within 
1 kpc of the sun (where the dust dominates). 

Using the iterative approach outlined in section 3, we 
looked for local minima in the reduced chi-square for a 
number of different local bubble sizes ranging from pc 
to 350 pc in steps of 25 pc. We found a best fit local 
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Fig. 10. — Our derived Ry values as a function of the absolute 
value of the Galactic latitude, 6, in bins of 5°. It appears that 
there may be a correllation between Ry and Galactic latitude. 



TABLE 1 
Best-Fit Dust Parameters 



Parameter 


Value 


Central Density 


47 ± 117 X 10-26 g 


Scale Height 


119 ± 15 pc 


Scale Length 


1.35 ± 5.25 kpc 


Local Bubble Radius 


150 pc 


Local Bubble Dust Fraction 


40 ± 26% 



bubble with a radius of 150 pc containing no dust. We 
later included the dust fraction in the local bubble as a 
free parameter in our model and found the best fit dust 
fraction in a 150 pc radius was 0.40, i.e. the local bubble 
contains 40% of the dust predicted by our exponential 
model. Because the local bubbles shape is highly irreg- 
ular, this procedure is only meant to roughly constrain 
the local underdensity of dust. Accordingly, the formal 
l-<j error in this fit is 26%. We performed an F test on 
our results to test whether the 40% dust bubble is statis- 
tically more likely than the 0% dust bubble, and found a 
significance of 0%. This result indicates that the addition 
of a parameter to constrain the dust fraction inside the 
bubble, although resulting in a lower chi-squared value, 
has no statistical significance. Our model therefore finds 
that a 0% dust bubble is as likely as a 40% dust bubble. 

The fitting procedure returned a dust scale height of 
152±33 pc. However, the small amount of data at low 
latitudes and our treatment of the dust structure in- 
side our spherical local bubble limits the precision of our 
method. As an alternative, we removed all stars with 
Ay = from the dataset before binning and fitting for 
the dust scale height, because we would not expect lines 
of sight without dust to scale exponentially. This method 
effectively removed any effects of the zero extinction lines 
of sight from our derived scale height. From this proce- 
dure, we found a dust scale height of 119±15 pc. The 
actual scale height is likely bracketed by the results from 
our two methods. Our derived scale length was nearly 
infinite for the first method, but after removing all stars 
with Ay=0, we found a value of 1.35±5.25 kpc. The 
large uncertainty in this parameter is due to a lack of 
high-quality spectra at significantly different values of 
Galactocentric radii. 

Figures [TT] and [12] show the theoretical extinction and 
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Fig. 11. — The theoretical Ay (color coding) assuming an expo- 
nential dust model that scales both in R and Z with a spherical 
local bubble that also scales exponentially but with 40% of the 
dust density. We can see from this model compared to our actual 
data (Figure O that the dust is considerably more patchy than a 
uniform exponential scaling. 
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Fig. 12. — The dust densities from our theoretical model (color 
coding), which finds a dust scale height of 119 ± 15 pc, a scale 
length of 1.35±5.25 kpc, and a spherical bubble with a radius of 150 
pc containing 40% of the dust predicted by a smooth exponential 
model. 

dust distributions, respectively, using the scale length 
and height from the second method above and the local 
bubble size and dust fraction from the first method. Ta- 
ble [1] contains the list of best-fit dust parameters used in 
these models. The actual extinction does vary consider- 
ably compared to the uniformly scaling best-fit model, 
which indicates the patchiness of the dust distribution. 
Figure [13] shows Ay divided by distance for the stars in 
our sample with Ay greater than 0. Ay divided by dis- 
tance is an indicator of dust because dust is proportional 
to the differential extinction. This means that extinc- 
tion divided by the distance is directly proportional to 
the average dust density along the entire line of sight to 
the star (see Equation 4). Figure [13] shows that the dust 
does appear to scale exponentially, on average, and fits 
our derived exponential extinction model very closely. 

4.4. Comparison with SFD 

Figure [14] shows a map in Galactic latitude and longi- 
tude of SFD's total column extinction subtracted from 
our Ry = 3.1 Ay values. Although our values were 
usually in close agreement to those of SFD, our values 
are often lower than the SFD values by a significant 
amount at lower Galactic latitudes. This is expected 
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Fig. 13. — Median Ay divided by distance, an indicator of dust, 
as a function of distance from the plane in bins of 25 pc. We have 
left out the values of Ay = 0, because we don't expect lines of 
sight with no detectable dust to scale exponentially. However, this 
plot shows that the lines of sight with significant dust tend to scale 
exponentially. Our best fit extinction model is plotted in red and 
closely matches the data. Note that the best fit extinction model 
is a fit to our extinction map and is not a fit to the data shown. 

due to the fact that we are not seeing through all of the 
dust near the Galactic plane, and our results fo r these 
regions are consistent with t hose of ISchlafiv et al.l (|2Q1Q| ) 
and S chlafiv fc Finkbeiner (2010). Areas where our re- 
sults are lower than those of SFD at higher latitudes 
are also likely due to stars that are in front of a sig- 
nificant fraction of the dust. Areas at higher latitudes 
that report higher Ay values than those of SFD values 
could be due to error in our extinction values, our as- 
sumption of Ry =3.1^ or possibly due to steep extinction 
gradients, where SFD m ay underestimate the extinction 
(|Arce fc Goodman] 1 1999[ ). Representative of this. Figure 
[T4] also shows comparisons with SFD for stars at dis- 
tances less than and greater than 200 pc. These maps 
show that our extinction values are more likely to be 
closer to those of SFD at greater distances. In the upper 
left of the middle map (distance < 200 pc) and in a few 
other lines of sight, however, we find values significantly 
higher than those of SFD. We will examine these features 
in more detail in future studies. 

5. CONCLUSIONS 

By fitting over 56,000 SDSS M dwarf spectra to the 
lFitzpatric]J(jl999[ ) extinction curve, we have created ex- 
tinction and Ry maps for the local Galaxy. Our results 
are summarized below. 

1. We developed a technique for empirically measur- 
ing the extinction and distance to stars for which 
we have SDSS spectroscopic data. Most previous 
techniques have used a two-color analysis to deter- 
mine extinction; our method uses the flux values 
across the spectral range from 5700 to 9200 A to 
determine the best possible extinction and distance 
fits, and so has the potential to be much more ac- 
curate. 

2. Using our fitting procedure, we deriv ed new dis- 
tances for most of the stars in the I West et al.l 
(f2011 ) SDSS DR7 M dwarf catalog, accounting for 
the proper amount of foreground dust extinction. 

3. We created an Ay map extending to approximately 
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Fig. 14. — A Galactic latitude and longitude map of the values of Ay determined by SFD subtracted from our values of Ay (color coding) 
for all the stars in our sample (left), stars with distance < 200 pc (center), and stars with distance > 200 pc (right). In this plot we assume 
an Ry =3.1 extinction law. Many lines of sight have significantly lower extinction than the SFD values, suggesting that these stars are in 
front of a significant portion of the dust column. We expect that the places where our values are lower than those of SFD at lo w latitudes 
are du e to the fact that we are not seeing through all of the dust near the plane. This is consistent with the findings of iSchlaflv et al.l 
(|2010l ) and Schlafiv & Finkbeiner (2010). In the center map (distance < 200 pc) there are a few more lines of sight where our values are 
lower than those of SFD, suggesting that these stars are in front of a significant portion of the dust column. There is an interesting feature 
in the upper left, where our values are higher than those of SFD at a variety of distances, which will be examined in more detail in future 
studies. In the map on the right (distance > 200 pc), there are few lines of sight where we find lower values than those of SFD, suggesting 
that many of these stars are beyond one dust scale height. 



2 kpc from the sun, generally at high Galactic lat- 
itudes. In the northern hemisphere, where most of 
the stars in the SDSS sample lie, the 6>-compressed 
extinction is relatively uniform. In the southern 
hemisphere, however, there are fewer SDSS lines 
of sight and so significant variations are apparent, 
showing that the three-dimensional dust distribu- 
tion often does not scale uniformly. 

4. Our extinction map compares favorably to the SFD 
map over the vast majority of the SDSS survey 
area. Our extinction values are often significantly 
less than those of SFD, most likely because we are 
seeing through only part of the total dust column. 
The median error in our extinction values is 0.029 
magnitudes. 

5. Within about 1 kpc of the sun, we deduced the 
value of Rv for many lines of sight, finding that 
the median Ry is 3.38 with a median uncertainty of 
0.42. The peak of our distribution of Ry (Figure[5]) 
is at a value of 3.1, in agreement with the generally 
accepted value for the diffuse interstellar medium. 
As we find no significant source of systematic error 
and a K-S test indicates that this distribution is 
inconsistent with a Gaussian of a = 0.42, we expect 
that the spread in our distribution indicates a real 
spread in Galactic Ry- 

6. We estimated the scale height of dust in the lo- 
cal Galaxy to be 119 ± 15 parsecs. This is less 
th an the value found by th e theoretical model 
of iDrimmel fc Spergell ()2001f ) which finds a value 
around 185 pc near the sun and roughly consistent 
with the HI model (long; th ougs:ht to be a tracer 
of dust) of Kalberla fc Kerp (2009) which finds an 
approximate value at all galactocentric radii of 150 
pc. We note that there is considerable variation in 
the structure of dust throughout the local Galaxy. 
In addition, we find some evidence for a spherical 
local bubble of radius 150 pc and with approxi- 
mately 40% of the normal dust density. We plan 



to better constrain the shape, size and dust con- 
tent of the local bubble through future work using 
nearby L dwarfs. The empirical data varies con- 
siderably compared our model, but on average it is 
in close agreement (Figure [T3|). Similarly, the true 
dust distribution is not nearly as uniform as our 
model. 

Our 3D map of Ay is most appropriate for estimat- 
ing extinction to stars within around 500 parsecs of the 
Galactic plane. Our map within this region, especially 
at small distances, has fairly high-resolution and low un- 
certainties, and using the SFD maps for these Galacto- 
centric heights could significantly overcorrect for extinc- 
tion. Due to an approximate (but not uniform) dust scale 
height of 119 parsecs, we can safely assume that stars 
further than 500 parsecs from the plane are behind es- 
sentially all of the dust column and that two-dimensional 
maps like SFD's higher resolution map are the most ap- 
propriate for these distances. 

In using our map of Ry, it is important to remember 
that our standard error tends to be around 0.42; however, 
regions of especially high and low Ry should be noted, 
especially for shorter wavelength observations where er- 
rors caused by assuming an average Ry can be signifi- 
cant. We expect that although for many lines of sight in 
the Galaxy Ry = 3.1 is a very good approximation, there 
is a broad spread in Galactic Ry. Because our Ry his- 
togram reports few Ry values near 2.1 or 5.5, we do not 
expect that Galactic Ry values are frequently outside of 
the range that we selected. However, when our selection 
effect is taken into account, the spread in Galactic Ry is 
even greater than that shown in our histogram. 

In the future, we plan to expand our extinction map 
by using the fitting procedure developed here for SDSS 
spectroscopic data from other stellar types. Our com- 
plete Ay, Ry, and distance dataset is accessible online in 
FITS format at http:/ /people, bu.edu/aawest /dust. html. 
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